Effects of the Fermi-surface shrinking on the optical sum rule in pnictides 
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In this paper we investigate the effects of the band shifts induced by the interband spin-fluctuation 
couphng on the optical sum rule in pnictides. We show that, despite the shrinking of the Fermi 
surfaces with respect to first-principle calculations, the charge-carrier concentration in each band is 
almost unchanged, with a substantial conservation of the total optical sum rule. However, a signifi- 
cant transfer of spectral weight occurs from low-energy coherent processes to incoherent ones, with 
practical consequences on the experimental estimate of the sum rule, that is carried out integrating 
the data up to a finite cut-off. This has profound consequences both on the absolute value of the sum 
rule and on its temperature dependence, that must be taken into account while discussing optical 
experiments in these systems. 



I. INTRODUCTION 

Since the discovery of superconductivity in iron-based 
superconductors a renewed interest emerged on the prop- 
erties of interacting multiband systems. Indeed, all the 
families of pnictides are semimetals, with several small 
hole and electron pockets at the Fermi level, originat- 
ing from almost empty hole and electron bands. Such a 
topology has been predicted by Density Funcional the- 
ory (DPT) calculations and confirmed by several exper- 
iments sensitive to the Fermi-surface structure, as de 
Haas-van Alphen^"'^ and angle-resolved photoemissioii 
spectroscopy (ARPES)'*"^. Despite the qualitative agree- 
ment with DPT predictions as far as the number and the 
character of the bands is concerned, the experimental 
data suggest the existence in the real materials of a band 
narrowing operating over an energy scale of hundredths 
meV and of a shrinking of the Permi-surface areas with 
respect to DPT. 

The first feature has been associated to the presence 
of electronic correlations^'^ leading to an enhancement of 
the band mass m of the carriers with respect to the DPT 
value TTiDFT, in agreement with several DFT+DMFT 
calculations giving m ~ 2mDFT- This effect can be 
possibly detected also by an experimental estimate of 
the optical sum rule, which relates the integrated op- 
tical conductivity iy{uj) to the carrier density n and 
band mass, W = duja{uj) oc n/m.^"'^^ Indeed, re- 
cent experiments in LaPePO show that the experimen- 
tal weight Woxp is significantly reduced with respect to 
the DPT value, Woxp/Wdft « 0.5.^^ Assuming that 
DPT correctly estimates the number of carriers n in 
each band, this result can be interpreted thus as an ef- 
fect of the band-mass renormalization mentioned above, 

Wexp/M^DFT muFT/m. 

In addition to the electronic correlations, operating 
on large energy scales, a retarded low-energy interac- 
tion, mediated by the spin-fiuctuations, is also present 
in pnictides. Within this context, the shrinking of 
the Permi-surface areas has been explained as an effect 
of the band shifts originating from the interband coupling 



of electron and hole bands mediated by low-energy spin- 
fluctuations.^^ It is worth noting that these Permi-surface 
shrinkings do not violate the total charge conservation, 
which, in multiband systems, is determined by the bal- 
ance of the charge carriers in all the bands. In contrast 
to the many-body picture, an alternative explanation for 
the Permi surface shrinkings has been proposed to be an 
intrinsic inaccuracy of the DPT calculations in determin- 
ing the top/bottom energy of the bands with respect to 
the Permi level. In this case, a simple rigid shift of the 
DPT bands could account for the observed shrinking of 
the Permi-surface areas^ without involving many-body 
interaction effects. A natural consequence of this shift 
would be a reduction of the carrier concentration in real 
materials with respect to the DPT prediction, with di- 
rect implications for the optical sum rule. Indeed, the 
experimental observation of a low value of Wcxp could 
not be attributed only to mass-renormalization effects, 
but should be associated (at least in part) to a reduction 
of n due to the rigid-band shift. ^^'^^ 

Motivated by this framework, in this paper we present 
an extensive analysis of the effect of the Permi-surface 
shrinking induced by the interband interactions on the 
charge-carrier conservation and on the optical properties 
of a multiband system. In particular we show that within 
this context, in contrast to a Permi-surface shrinking in- 
duced by rigid-band shift, the charge-carrier concen- 
tration in each band does not scale with the Permi area, 
and it is instead only weakly affected by its shrinking, 
leading to a redistribution of spectral weight between 
coherent and incoherent processes and to a substantial 
conservation of the optical sum rule with respect to the 
non-interacting case. Notably, as we will show, the redis- 
tribution of optical spectral weight is remarkably differ- 
ent with respect to the rigid band shift suggested in Ref. 
[16], where the reduction of the Permi area is refiected in 
a corresponding reduction of both the Drude peak and the 
incoherent band (Pig. 1). On the contrary, in the case of 
a Permi-surface shrinking induced by the retarded inter- 
band interactions, the reduction of the coherent Drude 
part associated to the Permi-surface shrinking is accom- 
panied by an additional transfer of spectral weight to 
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FIG. 1: (color online) Sketch of the optical conductivity in 
three representative cases: (a) infinite-band system interact- 
ing with a bosonic mode. The optical conductivity is charac- 
terized by a coherent Drude peak and an incoherent band. No 
Fermi-surface shrinking is present; (b) same as in (a), but with 
a Fermi-surface shrinking due to a rigid band shift. In this 
case both the Drude peak and the weight of the incoherent 
processes are reduced; (c) effects of a Fermi-surface shrinking 
due to a many-body interaction in a multiband system. In 
this case the total spectral weight is (almost) unchanged, but 
the Fermi-surface shrinking reflects in a reduction of the co- 
herent spectral weight in the Drude part, with an additional 
increase of the incoherent processes. 



incoherent processes (Fig. Ic). As we shall see, this has 
practical consequences on the experimental determina- 
tion of the optical sum rule, which is usually carried out 
integrating the data up to a finite cut-off. 

The structure of the paper is the following. In Sec. II 
we introduce the model and we show how in a multiband 
system with interband interactions the charge-carrier 
concentration in each band does not scale with the cor- 
responding Fermi-surface area. In Sec. Ill we compute 
the optical conductivity in the presence of interactions 
and we analyze the effects of the Fermi-surface shrinking 
on the optical spectra. In Sec. IV we focus on the tem- 
perature effects on the optical-conductivity sum rule, by 
considering a bosonic spectrum that is either constant in 
temperature or temperature dependent. In Sec. V we 
comment on the outcomes of our results for the exper- 
iments in pnictides, and in Sec. VI we summarize our 
conclusions. 



II. FERMI-SURFACE SHRINKING AND 
NUMBER OF CARRIERS 

In this paper we investigate the effect of a retarded 
interaction on the optical properties of a multiband sys- 
tem characterized by two peculiar features: (i) a strong 
particle-hole asymmetry of the bands and (ii) a predom- 
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FIG. 2: (color online) Sketch of the band structure and of the 
Fermi surfaces with the main interband interactions. 



inant interband character of the interaction. The out- 
comes of these features on the one-particle properties 
have been discussed in Refs. [15,18], and the main re- 
sults relevant to the present analysis will be recalled be- 
low. We will then follow here the same effective low- 
energy multiband scheme of Refs. [15,18], which contains 
also the main ingredients needed to discuss the physics 
of iron-based superconductors. In particular, we will fo- 
cus on the parameter values appropriate for LaFePO, 
where both de Haas-van Alphen measurements of the 
Fermi-surface shrinking^ and measurements of the op- 
tical conductivity^^ are available in the literature. 

We consider two hole-like bands (1,2) located close 
to the r point and two electron-like degenerate bands 
(3,4) close to the M point (see Fig. 2). All the bands 
close to the Fermi level can be schematized as parabolic, 
respectively 



for the hole bands and 
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for the electron ones. We assume also that carriers in- 
teract via a bosonic mode with a typical energy scale luq. 
In the case of pnictides, where the largest source of in- 
teractions comes from the exchange of spin fluctuations 
between quasi-nested hole and electron pockets, the inter- 
actions Vap between two bands a and /3 will have a pre- 
dominant interband character. As it has been discussed 
in Ref. ]18], the coupling to a low-energy mode cannot 
account for the overall reduction of the bandwidth with 
respect to DFT observed in ARPES experiments. This 
effect must be attributed to Coulomb repulsion occurring 
at high-energy scales,®'^ and we model it by assuming our 
input masses as a factor twice larger than the DFT esti- 
mate, nia ~ 2mQ^DFT-^^'^* To be more explicit, setting 
for convention the Fermi level = 0, we first estimate 
from DFT calculations the Fermi vectors fcp.a and the 
nearest band edge for each band, namely i^max,!, £'max.2 
for the hole-like bands, i?min,3) £'min,4 for the electron 
ones. From these we estimate the non-interacting mass 



for each band 
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sponding bare density of states = 7x1^0^ I^-kY? [a is 
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(eV) 


(cV) 


(cV-i) 


1 


1.16 


0.102 


-2.516 


0.382 


2 


2.28 


0.102 


-1.231 


0.750 


3,4 


1.58 


1.776 


-0.147 


0.520 



TABLE I: Microscopic band parameters used in this work. 
Following the approach of Ref. 15, we use DFT band parame- 
ters renormalized by a factor 2, to account for the correlation 
effects. 



here the in-plane lattice constant). Accordingly, the ef- 
fective band edge far from the Fermi level^^ follows from 
the relation A^^ = l/(-Emax,a — ^^min,a)- Finally, to ac- 
count for the correlation effects wc simply divide by 2 the 
values of the band edges extracted by DFT, which cor- 
responds to the mentioned increase of a factor 2 for the 
effective mass TOq and for the density of states N^^ . The 
estimated band parameters in the presence of correlation 
are summarized in Table I. 

In the presence of a retarded boson-mediated interac- 
tion, the Green's function Ga{z) for a generic band a can 
be written as: 



Ga(k,w) = 



1 



(3) 



where we introduced the short-hand notations 

Xa(w) = ReE„(w), r„(w) = -ImSo,(a;). 

The self-energy in the Matsubara space can be 

computed as 

S„(iw„) = -T^Va,i3D{0Jn-l^m)G^{iuJm), (4) 

where Gfi{ioJm) = J {(f\i/4:TT'^)Gi3{k,iojm) is the 
local one-particle Green's function, D{u)i) = 
J dQ2nB{fl)/{n'^ + bjf) is the propagator of the 
bosonic mode and B{fl) is the density of states of the 
bosonic excitations. Here we assumed for simplicity 
that the interaction does not depend on the momentum, 
so that also the self-energy is momentum-independent. 
Finally, the self-energy and the Green's function 

Gc(k, w) on the real- frequency axis can be obtained 
employing the standard Marsiglio-Schossmann-Carbotte 
analytical continuation^" and the number of charge 
carriers per band can be obtained as 



a = 1,2 



(5) 



for the hole bands and 



dujf{uj)N^{uj) 



a 



3,4 



(6) 



for the electron bands, where f{x) = l/[cxp{x/T) + 1] 
is the Fermi function, Na{uj) is the interacting density 
of states (DOS), iV(w) = iV° / (ieA(e, w), and where 
A(ek,w) = — (l/7r)ImG(k, w -|- iO+) is the one-particle 
spectral function. 

To account for a spin-mediated interaction mechanism 
we use here the Lorentzian spectrum^ ^ 



1 flUJQ 



(7) 



with a characteristic energy scale luq = 20 meV. From the 
above relations we can introduce the matrix of the dimen- 
sionless coupling X^^^ = Va,fiNpD{0), which is related to 
the low-energy mass renormalization due to the retarded 
interaction as m* = (1 + Xa)ma, where Xa = '}Zj3 ^a,i3- 
Following Ref. 15 we take V^^ = V for the (a, /3) values 
shown in Fig. 2, and Va^ = otherwise. This choice 
is appropriate for LaFePO, while for 122 pnictides the 
anisotropy of Va0 must be taken into account.^* Nonethe- 
less, all the results we will discuss below depend only 
quantitatively on the exact form of the spectrum and 
of the Va/3 matrix, and will be qualitatively the same 
also for other choices of the boson mediator and of the 
interactions, provided that interband terms Va^,a ^ (3 
dominate over the intraband ones. 

As we discussed in Ref. 15, the strong particle- hole 
asymmetry of each band in pnictides forces us to calcu- 
late the Eliashberg self energy taking into account the 
finite bandwidth. In contrast to the usual infinite-band 
Eliashberg approximation, this leads to a finite value of 
Xa(w = 0), and to a corresponding change of size of 
the Fermi surface which depends on the interband vs. 
intraband character of the interaction. In the case of 
dominant interband scattering, in particular, these finite- 
bandwidth self-energy effects result in a shrinking of the 
Fermi surface, whose new Fermi vectors can be obtained 
from the poles of the Green's function (3) as 
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3, 4) we have 



1,2), while for electron bands (a 
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Note that, as shown in Ref. 15, for interband interactions 

Xa (0) is negative for the hole bands while it is positive for 
the electronic one, so that k-p^a is reduced in both cases. 
From the Fermi area we can define a coherent charge- 
carrier concentration per band: 



Ua = 



27r ' 



(10) 



where we have taken into account the spin degeneracy 
Ng = 2. In the case of a rigid band shift, as hypothesized 
in Ref. 16, is also equivalent to the total charge- 
carrier concentration per band n^. Things are however 
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FIG. 3: (color online) Variation of the hole and electron den- 
sities as a function of the energy shift at the Fermi surface. 
Here Ua is the total number of carriers in each band, while ha 
is the number of coherent carriers related to the Fermi area 
according Eq. (10) 



deeply different when the shrinking of the Fermi surface 
is induced by a retarded interaction mediated by a low- 
energy boson. Indeed, in this case although the reduction 
of the Fermi areas and hence of ha is remarkable, the to- 
tal charge carrier concentration per band is almost 
unaffected. Note that this result does not violate the 
Luttinger theorem, since for a multiband system there 
is no equivalence within each band between the charge 
concentration n„ and its corresponding Fermi area en- 
coded in ha, as it is instead the case in a single-band 
system. Indeed, the usual argument used to prove such 
an equivalence for the single-band case^'^ shows that in 
the multiband case the Luttinger theorem only applies to 
the sum of the carrier density in all the bands. In Fig. 3 
we plot the variation of both the total Ua and the coher- 
ent ha carrier density for each band as a function of the 
corresponding band shift Xq(0), which is proportional to 
the coupling to the bosonic mode.-"^^ The corresponding 
shrinking of the Fermi area can be directly obtained by 
Eqs. (8)-(9). As one can see, a band shift Xa{0) ~ 25 
meV (that corresponds to the Fermi-surface shrinking ob- 
served experimentally in LaFePO^'^^) leads to a remark- 
able reduction of ha, up to almost 25%, whereas the 
reduction of Ua is much less pronounced, with largest 
variations of Ua of the order of 5%. 

In order to understand why the total charge-carrier 
concentration per band is almost unaffected by the 
interaction in spite of the strong shrinking of the Fermi 
areas, let us examine the behavior of the interacting DOS 
N(uj). As an explicative example, we compare in Fig. 4 
the DOS of the electron band 3 in three cases: (a) with- 
out interaction, (b) in the presence of a rigid-band shift 
equivalent to X3(0)i ^'^'^ (c) in the presence of the cou- 
pling to spin fluctuations. Here we used a relatively large 
value of the coupling {V — 0.92 eV, leading to X3(0) ~ 60 
meV) to make more visible in the figure the effects of the 
interaction. In each panel the shaded area corresponds 
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FIG. 4: (color online) DOS of the electron band 3 in three 
representative cases: (a) bare DOS; (b) rigid-band shift; (c) 
interacting case, computed using V — 0.92 eV. The shaded 
areas correspond to the particle number at T = in the three 
cases. The arrow in panel (c) marks the new band edge due 
to the energy shift induced by X3(0)- Note however that the 
DOS extends also below this limit due to the finite value Fa (u) 
of the imaginary part of the self-energy. 

to the integral of the DOS up to Fermi level (i.e. uj — 
in our notation) , which gives the carrier density at T = 
according to Eq. (6). In the case of a rigid-band shift the 
number of carriers decreases following the Fermi-surface 
shrinking (10). However, as shown in Fig. 4c, when the 
shrinking is due to a many-body interaction the redis- 
tribution of spectral weight in the DOS is remarkably 
different. Here the presence of a finite Xsi'^) is reflected, 
in analogy with the rigid-band case, in a shift of the band 
bottom with respect to the non-interacting case, visible 
at the energy level where the quasi-particle DOS has a 
rapid drop (blue arrow in Fig. 4c). However, this effect 
is balanced by a strong redistribution of spectral weight 
in an extended tail of the DOS for uj < ii^min.s = —0.15 
eV, which (almost) compensates the band-edge shift and 
explains why is only slightly smaller than the bare 
value, as shown Fig. 3. Such long tail in the DOS is due 
to the incoherent states induced by the imaginary part of 
the self-energy r3(a;), and it is a characteristic signature 
of the interband nature of the interaction. 

To better understand this issue, we plot in Fig. 5 both 
the DOS and the self-energy for the band 2 and 3 over a 
larger energy scale, which allows us to show also that sim- 
ilar effects occur for both electron- and hole-like bands. 
In panel Fig. 5b we compare the corresponding imagi- 
nary parts of the self-energy Ta{uj). It is important to 
notice here that in the hole- like band r2(a;) has a finite 
support that coincides approximately with the DOS of 
the electronic band N3{uj). The same occurs for r3(LLi), 
which is controlled by the support of iV2(w). This is 
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FIG. 5: (color online). Panel (a): DOS of the hole band 2 
and of the electron band 3 (solid lines), displayed along with 
the bare DOS (dashed lines). To better resolve the two bands 
we used a negative y axis for the DOS of the band 2. The 
shaded areas give the number of carriers Ua in each band at 
T = 0, according to Eqs. (5)-(6). Panels (b) and (c): energy 
dependence of the imaginary and of the real part of the self- 
energy for the same bands. Note that V2{ijj) is different from 
zero in an energy range corresponding to the DOS of the band 
3, and viceversa. 



a direct consequence of the interband character of the 
interaction, where the energy support of [to) is deter- 
mined by the DOS of the band (3 to which the band a 
is coupled. This situation has drastic effects on the re- 
sulting DOS. Indeed, as shown in Fig. 5a, the incoherent 
states induced by Ta{uj) are thus created mainly outside 
the energy range of the band a, namely below the band 
bottom for the electron dispersions and helow the band 
top for the hole ones, and they are responsible for the 
long tails in the corresponding DOS. It is worth stress- 
ing here that such tails are not due to a strong-coupling 
effect, but simply to a predominant interband nature of 
the interaction. 



III. OPTICAL CONDUCTIVITY AND 
SPECTRAL WEIGHT 

In the previous Section we have seen how the Fermi- 
surface shrinking generated by a retarded interaction is 
not directly reflected in a substantial reduction of the 



charge carrier density Ua for the corresponding band, al- 
though the total number of coherent states ha is strongly 
affected by that shrinking. As we shall see in the present 
Section, such different behavior can be properly detected 
in the analysis of the optical properties. In particular we 
show that the optical sum rule, which is directly related 
to total number of carrier Ua in each band, is almost 
unaffected by the Fermi-surface shrinking, whereas the 
DC conductivity and the spectral weight in the coherent 
Drude peak scale with ha, so that they are strongly re- 
duced by the Fermi-surface shrinking, with a transfer of 
spectral weight to the intraband incoherent processes. 

To address this issue we compute explicitly the op- 
tical conductivity in the presence of the retarded in- 
terband interaction. Since for a k-independent self- 
energy vertex corrections vanish, (j{uj) can be computed 
in the simple-bubble approximation. We can write thus 
(7{uj) = Y,a o'a(^) ' where 



a a (w) 



2^ ^^/(z-M + c.)-/(^-/i) 

h 



where Wk,a = {X/K)d£-k_^a/ dkx is the quasiparticle velocity 
and the spectral function can be written explicitly as: 



7r[c.-£-XaM]2 + [rr(c.)] 



(12) 



Here we considered also a finite contribution of disor- 
der to the total quasi-particle scattering rate, r^P(a;) = 
r„(w)-Hro. Since as T-J>0r„(0) ^ Owe choose Fq = 10 
meV by estimating the approximate width of the low- 
energy optical spectra from Ref.^^ (see also Eq. (20)-(21) 
below). By converting the k integration in an energy 
integration and considering the parabolic hole-like and 
electronic-like dispersions described in Eqs. (l)-(2), we 
obtain for the hole bands (a = 1, 2) 



(7 a (t^) 



dz 



f{z - fi + uj) - f{z - ^) 



X / de{E^s^^,a-e)Aa{s,z + uj)Aa{e,z) 



(13) 



while for the electron bands [a ~ 3,4) 



cra(w) = — - / d: 

n- J-oo 



J{z - fl + Oj) - f{z - fl) 



UJ 



X / dt{e ~ E^i^^a)Aa{£,z + uj)Aa{e,z). 

(14) 

The resulting optical conductivity (t(w) at low temper- 
ature is shown in Fig. 6 for the same coupling value 
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FIG. 6: (color online) Left axis: Frequency dependence of 
the optical conductivity at T = 5 K as obtained from (13)- 
(14) (solid line), and low-energy Drude-like peak as obtained 
from Eq. (21) (dashed line). The small discrepancy for the 
dc values stems from the use of the approximate formulas 
(19)-(18) in the Drude term. To permit a simpler comparison 
with the experiments we divided the two-dimensional results 
by the interlayer distance d — 8.5 A. Right axis: integrated 
spectral weight W(lS), compared to the integral of the Drude 
part only, W-q^lo). The arrows mark the bare value Wo, the 
asymptotic value W (15) in the interacting case, and the value 
W corresponding to the reduced Fermi surfaces. 



To enlighten this issue and its relation to the Fermi- 
surface shrinking let us first evaluate the dc conductiv- 
ity in the zero temperature limit, where we can write 
in Eqs. (13)-(14) lim„^o[/(^-M + ^)-/(^-M)]/w « 
— (5(z — /i). By performing the remaining integration over 
e analytically we get thus: 



arctan 
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(16) 
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(17) 



for the holes. For T — > the only source of damping 
is due to disorder, so that r''P(0) reduces to impurity 
scattering Fq. Since Fg is usually much smaller than the 
distance of the both/top of the band from the Fermi level 
[including also the shift x(0)], the above expressions sim- 
plify considerably. Indeed, using the expressions (8), (9), 
(10), we obtain for the hole bands (a = 1,2) 



V — 0.46 eV used in Ref. 15 to reproduce the Fermi- 
surface shrinking measured by de Haas-van Alphen ex- 
periments in LaFePO. Also shown in the same plot is the 
integrated optical spectral weight W{u!) — /g" duj'a{uj') 
as a function of the cut-off frequency. When the integra- 
tion frequency goes to infinity the total spectral weight 
is given by the optical sum rule, that for a parabolic- 
band approximation as the one we are using here simply 
reduces to^*''-^^ 



W= I doja{Lj) = ^-Y,^- (15) 



This asymptotic value is marked in Fig. 6 by a solid ar- 
row, along with the sum-rule value Wq in the absence 
of the coupling to spin fluctuations. As we can see, the 
discrepancy between W and Wq results to be of order 
of 5 %, pointing out that the effect of the Fermi-surface 
shrinking induced by the many-body interaction is quite 
small in the total sum rule. However, by direct inspec- 
tion of Fig. 6 one can also see that even at frequencies 
w ~ 3000 cm^^, corresponding to about 20ujq, the col- 
lected spectral weight W{ll!) still differs by about 20% 
from its asymptotic limit (15), and it is near instead to 
the spectral weight corresponding to a rigid-band shift, 
i.e. W — Tre^ na/2ma- Such a lack of saturation of 
W{lu) over frequency scales much larger than the typical 
boson energy luq is unusual for single-band systems, ^^'^'^ 
and indicates that for the multiband case discussed here a 
large spectral contribution is allocated at relatively high 
energies. 



O'dc 



Xo(0)] = 



27r?iFqp(0) ' 
and for the electron bands [a — 3,4) 



27r?iFqp(0) 



Xa(0)] - 



(18) 



(19) 



where we introduced the effective mass to* = (1 + Aq)toc 
and the transport scattering time 



2rr(o) 



(20) 



As it is evident from Eqs. (18)-(19), even though the 
total carrier numbers are almost unaffected by the 
shrinking of the Fermi areas (see Fig. 3), the dc conduc- 
tivity is instead controlled by the coherent carrier density 
fia ■ More generally, one can generalize this result to finite 
frequency to get an approximate expression of the coher- 
ent Drude peak a^{uj) in the presence of the many-body 
interaction as 



1 



to: 



1 + {ojTtr.a) 



(21) 



The coherent Drude peak described in Eq. (21) is also 
shown in Fig. 6, along with its integrated spectral weight 
Wd on the right-side scale. The comparison with W{lu) 
shows that Wd(ll') saturates quite rapidly to its asymp- 
totic value Tre'^n/2'm* , while W{uj) steadily increases in 
this range of energies, showing that all the remaining 
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spectral weight needed to recover the value (15) must 
originate from incoherent processes at higher frequen- 
cies. Note that this result has to be contrasted with 
the standard infinite-band case, where the real part of 
the self-energy vanishes and the spectral- weight distribu- 
tion is only due to the mass renormalization, so that the 
coherent Drude part Wd oc ric/m* (see Fig. la).^^ On 
the contrary, here the many-body interaction gives rise to 
an additional Fermi-surface shrinking which is refiected 
in a further reduction of the coherent spectral weight 
W oc fia/m'^ (see Fig. Ic). 



IV. TEMPERATURE DEPENDENCE 

In the previous section we have seen that the many- 
body interband interaction plays a non trivial role on the 
spectral properties of systems with strongly particle-hole 
asymmetric bands, and in particular in pnictides. This 
is pointed out for instance in the coherent contribution 
to the optical conductivity described in Eq. (21). In infi- 
nite bandwidth systems, assuming the bosonic spectrum 
B{Q) to be weakly temperature dependent, the main 
temperature dependence of the optical conductivity in 
Eq. (21) is through the parameter Ttr.a{T) [or equiva- 
lently r^P(T)] whereas h — n is temperature indepen- 
dent. This scenario has to be contrasted with the case 
of a multiband system with strong particle-hole asymme- 
try, where the real part of the self-energy Xa{0), which 
gives rise to the Fermi surface shrinking, is itself temper- 
ature dependent and it is reflected in a coherent charge 
concentration fia which acquires a finite temperature de- 
pendence. Such will be analyzed in the Section 
IV A, whereas in Section IV B we consider the addi- 
tional effects due to a temperature-dependent bosonic 
spectrum, as suggested by a number of experiments in 
pnictides. In this case the dimensionless electron-boson 
coupling Aq and the low-energy mass renormalization 
^a/n^a = 1 + acquirc their own temperature depen- 
dence, which can compete with the previously discussed 
effects. As we shall see, transport and optical properties 
can present quite different temperature behaviors in the 
two cases. 



A. Temperature-independent bosonic spectrum 

We consider first the case where the bosonic spectrum 
is not itself temperature dependent. In this case the for- 
malism is the same as we discussed so far, and we should 
only account for the temperature evolution of the self- 
energy. In Fig. 7 we show as an example the temper- 
ature dependence of the real and imaginary part of the 
self-energy at zero frequency for band 2. As we can see, 
the absolute value of both quantities increases as the tem- 
perature increases, due to the thermal excitations of the 
bosonic mode. To get some analytical insight into these 
results one can employ the first-order expansion for the 
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FIG. 7: (color online) Left axis: temperature dependence of 
the real and imaginary part of the self-energy at zero fre- 
quency for band 2. Right axis: temperature dependence of the 
coherent carrier density n2 obtained from the Fermi-surface 
area according Eq. (10). 



self-energy, computed using the non-interacting Green's 
function in Eq. (4). The corresponding expression for 
Xq(0,T = 0) has been derived in Ref. 15 in the case of 
an Einstein bosonic spectrum, B{n) = {uj-e/2)S(CI — uj-e). 
By performing the same calculation at finite T one can 
see that the most relevant temperature corrections are 
related to the thermal excitation of the boson: 



xf(0,T) 



In 
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1 + ^bi.,/T) 

EcJ3 



(22) 



where Ec,p = min(£'inax,/3j ^'min,^) is the nearest band 
edge, and b{x) = l/[exp(a;/T) — 1] is the Bose function. 
In the case of a spin-fluctuation spectrum the direct com- 
parison with the results in Fig. 7 shows that a similar 
expression holds, provided the identification of the char- 
acteristic energy scale we ~ 2a;o- The imaginary part 
of the self-energy at zero frequency can be also com- 
puted exactly in the same approximation for tempera- 
tures T <^ Ec- 



rf(0,T) 



/3 



(23) 



An analytical results can be derived for Ta{0, T) also for 
the spin-fluctuation spectrum, and we obtain 



ra(o,r) 



-2rarctan — > \ap^ 



(24) 



which is in good agreement with the results of Fig. 7. 

For what concerns the optical spectrum, the ther- 
mal excitation of the bosonic mode has then two con- 
sequences: (i) a further reduction of the Fermi-surface 
areas, encoded in the decrease of na{T) with tempera- 
ture, as shown in Fig. 7; (ii) a progressive broadening of 
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FIG. 8: (color online) Optical conductivity for three represen- 
tative temperatures T = 10, 100, 200 K using a temperature- 
independent bosonic spectrum _B(f2) and an impurity scatter- 
ing rate To = 10 meV. Inset: temperature dependence of the 
integrated spectral weight W{ujc) with ujc = 3000 cm~^. 



the optical conductivity features, due to the increase of 
TaiT). While the latter effect is customary also in the 
single-band infinite-bandwidth case, the former one is pe- 
culiar of multiband systems with interband-dominating 
interactions, and leads to an unusual temperature de- 
pendence of the partial optical sum rule W{ujc)- Indeed, 
as we have seen in Section III, since the Fermi-surface 
shrinking gives rise to a transfer of spectral weight from 
the low-energy Drude peak oc n to the incoherent struc- 
ture extending up to high energies, a finite temperature 
dependence of h reflects directly in a sizable temperature- 
dependent depletion of W{ujc) when a finite cut-off Uc is 
employed. 

In Fig. 8 we show the evolution of the optical conduc- 
tivity for three different temperatures, using again the 
multiband scattering matrix Aq^ and the bosonic spec- 
trum estimated in Ref. 15 to reproduce the Fermi-surface 
shrinking of LaFePO at T w 10 K. Here we assume the 
all these parameters are temperature independent and 
we also use a constant impurity scattering rate Tq — 10 
meV. The corresponding variation of the spectral weight 
W{lUc) as a function of the temperature is reported in 
the inset, showing an overall reduction of about 10 %. 
Note that such a variation is an order of magnitude larger 
than the one expected in a single-band system for the 
same coupling value, even when the presence of a finite 
bandwidth^"'''^^ and a finite cut-off^^^^^ are taken into ac- 
count. Quite remarkably this temperature variation is 
as large as the one observed in cuprate superconductors, 
where it has been interpreted as a consequence of strong 
correlations. ^'''^^^^^ In our case however this effect comes 
primarily from the temperature dependence of the coher- 
ent carriers number, that is relatively large even when the 
coupling to the bosonic mode is not particularly strong. 



B. Temperature-dependent bosonic spectrum 

As we have seen in Sec. Ill, the low-energy part of the 
optical spectrum resembles a renormalized Drude model, 
with an inverse scattering time proportional to the 
imaginary part of the self-energy. As a consequence, if 
the boson spectrum is constant in temperature the broad- 
ening of the optical conductivity as the temperature in- 
creases is directly proportional to the increase of rQ(0), 
as shown for instance in the case of Fig. 8. With this 
choice, however, we get an optical width of the Drude- 
like low-frequency part at room temperature of about 600 
cm~^, much larger than then scattering rate t^.^ ~ 320 
cm~^ that can be fitted from the available experimen- 
tal data at T = 300 K.^^ This observation suggests thus 
that the increase of rQ(0) with temperature is smaller 
than predicted by Eq. (23). To account for this discrep- 
ancy we investigate here the possibility that the bosonic 
spectrum is itself temperature dependent, as it has been 
suggested also by a recent analysis of optical spectra in 
122 compounds. In particular, in Ref. 28 it has been 
shown that the bosonic spectrum extracted from optics 
hardens as T increases, while its total spectral weight, 
encoded in the dimensionless couplings A^^, decreases. 
This trend is also consistent with direct measurements 
of the spin-fiuctuations spectrum by means of neutron 
scattering in 122 compounds. Following Ref. 29, we 
model thus such behavior by assuming a temperature- 
dependent bosonic spectrum of the form: 



B{n) 



h{T) ujquj 



where 



and 



h{T) = 



T + Te' 



(25) 



(26) 



(27) 



We chose the value of Tq — 150 K in order to repro- 
duce the experimental value of the Drude peak in the 
optical data of LaFePO at T = 300 KP The evolution 
of the spectrum B{^1) for three representative tempera- 
tures T = 0, 100, 300 K and the corresponding temper- 
ature dependence of the coupling constant are shown in 
Fig. 9a and 9b, respectively. Focusing for simplicity 
once again on band 2, we plot in panel (c) the resulting 
real and imaginary parts of the self-energy X2(0)j ^2(0). 
As one can see, for this choice of parameters the reduc- 
tion of A2(T) is strong enough to result in a decrease of 
1X2(0)1 with the temperature and in a smaller increase 
of r2(0). As a consequence, also the temperature evo- 
lution of the optical conductivity, shown in Fig. 10, is 
markedly different than before. In particular, the reduc- 
tion of |xa(0)| with the temperature is reflected now in 
an increase of na{T) (see Fig. 9c). The spectral weight 
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FIG. 9: (color online) (a) Temperature evolution of the 
bosonic spectrum according to Eqs. (25)- (27). (b) Tempera- 
ture evolution of the dimensionless coupling constant in each 
band, (c) Corresponding temperature dependence of X2(0) 
and r2(0) (left axis) and of (right axis) in band 2 using 
the temperature-dependent bosonic spectrum i3(f2) of panel 
(a). 
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FIG. 10: (color online). Optical conductivity for three rep- 
resentative temperatures using the temperature-dependent 
bosonic spectrum _B(f2) shown in Fig. 9a. Here we used a 
constant impurity scattering rate Fq = 10 meV. Inset: tem- 
perature dependence of the integrated spectral weight W{loc) 
with ujc ~ 3000 cm~^. 



moves thus progressively toward the low-energy coher- 
ent part of the spectrum, leading also to an effective in- 
crease of the spectral weight W(yi^ integrated in a fixed 
frequency windows, in sharp contrast to what expected 
in ordinary single-band (interacting and non-interacting) 
systems. ^'^'^^'^^ 



V. COMPARISON WITH THE EXPERIMENTS 

Even though a detailed comparison with the experi- 
mental optical spectra is beyond the scope of the present 
paper, we would like to comment on the possible conse- 



quences of our findings for the physics of pnictides. A 
first issue concerns the role of the cut-off in the estimate 
of both the absolute value of the spectral weight and its 
temperature dependence. From the experimental point 
of view, a finite cut-off is always employed in the esti- 
mate of the sum rule to avoid inclusion of interband tran- 
sitions, so that Wcxp ~ 'WiuJc)- In pnictides the choice 
of Wc is dictated by the observed threshold to interband 
transitions, that occurs usually around uic = 2000 — 3000 
cm~^. This is indeed the case in Ref. [12], where the ex- 
perimental estimate Woxp for LaFePO includes only pro- 
cesses up to a cut-off frequency of — 3000 cm~^. As we 
have shown in Fig. 6, the calculated spectral weight col- 
lected up to this same cut-off is about 20% of its asymp- 
totic value W , which is instead almost unchanged (within 
^5%) from the bare value Wq (we remind here that in 
our model Wq. While comparing with the experiments 
we must recall that our starting model has renormalized 
DFT bands, so that Wq — O.SM^dft- As a consequence, 
after inclusion of the low-energy spin fluctuations we es- 
timate overall W{y)c) — 0.4Wdft, i-e. the main source 
of spectral-weight reduction still comes from correlation 
effects. We observe however that for the range of param- 
eters used here W{ujc) is very near to the spectral weight 
W that one would obtain from a rigid- band shift. The 
reason is that, as we have seen in Sec. II, the incoherent 
part of each spectral function extends up to energy scales 
of the order of the bandwidth of the coupled band, i.e. 
over energy scales much larger than the typical boson 
energy ujq. As a consequence, even though there is no 
theoretical sum rule relating W{uJc) to W, on the energy 
scales considered in the experiments this could be a good 
approximation for it. Note however that the situation 
could be slightly different in the 122 compounds, where 
the most populated band is also the less coupled to spin 
fluctuations,^^ making it difficult to predict a-priori in a 
precise quantitative way the relevance of the shrinking 
on the absolute value of optical sum rule with respect to 
DFT. 

A more direct effect of the Fermi-surface shrinking on 
the optical data is instead the temperature dependence 
of sum rule. Indeed, as we discussed in the previous Sec- 
tion, the temperature dependence of W(wc) is strongly 
affected by the temperature evolution of the coherent 
charge density fia, which in turn is an indirect probe of 
the temperature evolution of the spin-fluctuation spec- 
trum. Quite interestingly, recent optical measurements 
in both carrier-doped"^*^'"^^ and isovalent-substituted'^^ 122 
compounds show that Woxp increases with increasing 
temperature, in contrast to any ordinary behavior in 
single-band system. ^"'-"^^'^^ Within our scenario this sug- 
gests that the coupling to the spin fluctuations is actually 
temperature dependent, leading to an increase of TZfj as 
the temperature increases, see Fig. 9c, and to a progres- 
sive population of the low-energy optical spectrum, see 
Fig. 10. It is worth noting that our approach completely 
neglects the role of optical (q=0) interband transitions, 
that are presumably less relevant than the single-band 
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transitions in the range of frequencies considered here. 
However, a detailed analysis of the role of these inter- 
band transitions and of their possible relevance in the 
experiments is at the moment an open question that de- 
serves further theoretical investigation. 

A third outcome of our results is related to the possible 
observation of these anomalous temperature behaviors in 
the dc transport properties. Indeed, from the approxi- 
mate expressions (18)- (19) for the dc conductivity it is 
clear that as the temperature increases both ha{T) and 
Ta{T) contribute to the temperature dependence of the 
dc conductivity. In other words for a multiband system 
with dominant interband interactions one cannot simply 
attribute the T dependence of the resistivity to scatter- 
ing processes, since also the density of coherent carriers 
can have a non trivial temperature dependence. This ob- 
servation could shed new light also on the comparative 
analysis of dc and Hall conductivities carried out in Refs. 
[33,34], where part of the anomalous behavior reported 
there could be accounted for by relaxing the assumption 
of a constant density of carriers with the temperature. 

VI. CONCLUSIONS 

In this paper we analyzed systematically the behavior 
of the optical-conductivity sum rule in a multiband sys- 
tem with dominant interband interactions, as it is the 
case in pnictide superconductors. Due to the strong 
particle-hole asymmetry of these systems, the coupling 
to a collective mode leads to a many-body induced band 
shift and to a corresponding shrinking of the Fermi- 
surface areas. In this context, we have shown however 
that the number of particles in each band does not scale 
with the Fermi areas, but it is almost unchanged. Such 
behavior has been explained in terms of a transfer of the 
spectral weight from the coherent states related to the 
reduced Fermi areas to the unoccupied part of the non- 
interacting bands. As we have discussed in the present 
paper, the effects of such redistribution of spectral weight 



on the optical spectrum are quite different from what ex- 
pected in single-band systems. In particular, we have 
shown that only a fraction ~ h/m* of the total weight 
~ n/m is collected in the coherent low-energy part of the 
spectrum, while the rest is recovered as higher-energy in- 
coherent processes. As a consequence, the unavoidable 
use of a finite cut-off oJc in the experimental determina- 
tion of the optical sum rule is expected to have a much 
stronger effect than in the single-band case. There are 
two main consequences: (i) W^xp usually underestimates 
the theoretical sum-rule value W in a percentage that 
cannot be established a-priori, since it depends on the 
parameter values. In the case of LaFePO discussed here, 
Wexp ttirns out to be about 20% smaller than W. As a 
consequence, if one assumes that W is already strongly 
reduced (~ 50%) by correlations,^'^ the many-body ef- 
fect due to spin fluctuations is a small fraction (~ 10%) 
of the DFT value Wdft', (ii) the temperature evolu- 
tion of the optical sum rule turns out to be an indirect 
measurement of the temperature evolution of the Fermi- 
surface areas. For instance, a temperature dependence 
of the spin-fluctuation spectrum could suggest that the 
cj-integrated spectral weight W{ujc) could increase with 
temperature, in sharp contrast to the case of interacting 
single-band systems, as for example cuprate supercon- 
ductors. All these findings pose strong constraints on 
the analysis of optical data in Fe-based materials, show- 
ing once more that the peculiar semimetal and multiband 
character of these systems forces us to revise standard 
paradigms appropriate for single-band correlated mate- 
rials. 
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